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1 INTRODUCTION 



The 8 Set stars are low-overtone pressure mode pulsators 
on or near the main sequenc e between s pectral types A2V 
to FOV and A3III to F5III (|Kurt zll2000l ). whose pulsation 
frequencies range from 3 to 80 d -1 (i.e., 18 min to 8 h in pe- 
riod). Their pulsations are driven by the opacity (k) mech- 
anism operating on the Hell partial ionisation zone. These 
intermediate mass (1.5 — 2.5 Mq) stars are mostly Pop. I 
stars, but some Pop. II 8 Se t stars do exist and are known 
as SX Phe variables (see, e.g. iTempleton. Basu fc Demarquel 
120021 ) . 



The 5 Set instability strip is located at the junction be- 
tween the classical Cepheid instability strip and the main 
sequence. It is not populated solely with 8 Set variables with 
normal spectra, however. The <5 Set stars have as their neigh- 
bours: the photometrically less variable, chemically peculiar, 
classical metallic-lined (Am) stars; the marginal Am (Am:, 
pronounced "Am colon") stars with slightly milder abun- 
dance anomalies; the early or hot Am stars; the pPup stars 
(evolved Am stars); the Ap stars (where the 'p' means pe- 
culiar); the 7 Dor stars; the A Boo stars; and the 'normal' A- 
stars. This paper concerns KIC 3429637, a pPup star show- 
ing 5 Set pulsations. 



ABSTRACT 

KIC 3429637 (HD 178875) is a 8 Set star whose light-curve shows continuous pulsa- 
tional amplitude growth in Kepler Mission photometry. Analysis of the three largest 
amplitude peaks in the Fourier transform indicates different growth rates for all three. 
We have ruled out instrumental causes, and determine the amplitude growth to be 
intrinsic to the star. We calculate time-dependent convection models and compare 
them with the observations. We confirm earlier characterisations that KIC 3429637 is 
a marginal Am star through the analysis of new spectroscopic data. With the data 
presently available, a plausible cause of the amplitude growth is increasing pulsational 
driving as evolutionary changes shift the Hell driving zone deeper in this pPuppis star. 
If this model is correct, then we are watching real-time stellar evolutionary changes. 

1.1 The Am stars 

The classical Am stars are those whose Call K-line types 
appear too early for their hydrogen line types, and metallic- 
line types appear too late, such that the spectral types in- 
ferred from the Call K- and metal-lines differ by five or more 
spectral subclasses. Such stars are therefore commonly given 
three spectral types, e.g. AmkA8hA9mF3, corresponding to 
the K-line, hydrogen-lines and metallic-lines, respectively. 
The Am group makes up a si gnificant fraction of late A 
stars - up to 50 per cent at A8 (|Smith|[l973l ; ISmallev et alj 
l201ll ) , a nd feature s abundance anomalies typically of factors 
of ±10 (Abt 200£j). Historically, the classical Am stars were 
thought not to pulsate, but with the micromagnitude pho- 
tometric precision of the space-based Kepler mission six of 
the ten known A m stars in the Kepler field of view are ob- 
served to pulsate |Balona et alj l2011al ) , and ISmallev et alj 
(|201ll ) found 200 of 1600 metallic-lined A and F dwarfs to 
be pulsating from SuperWASP data. 

The Am stars rotate sufficiently slowly that turbu- 
lent motions do not prevent gravitational settling of helium 
jBaglin et al.lfl973l ). When enough helium is drained from 
the Hell partial ionisation zone, the star no longer pulsates. 
Meanwhile, ions with absorption lines near the peak wave- 
length of the photon energy distribution are radiatively lev- 
itated towards the surface, and others, notably He, C, Ca 
and Sc, gr avitationally settle, thus providing the abundance 
anomalies (|Kurtzll 19891 ). 
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The marginal Am stars show fewer than five spectral 
subclasses between the Ca II K and metallic lines. Such stars 
are therefore less extreme example s of the classic al Am stars, 
with milder abundance anomalies. iKurtj (|l97Sf ) showed the 
majority of the Am: stars lie at the blue edge of the 8 Set 
instability strip (his Fig. 2), which makes them difficult to 
distinguish from the 'hot Am stars' whose primary charac- 
teristics are Am anomalies in stars hotter than A4. Marginal 
Am stars are also found at the cool edge of the instability 
strip, however, with spectra similar to the his t orical <5Del 
stars. We refer the reader to iGrav fe Garrisonl jl989t) for a 
discussion of why the <5Del classification has been dropped 
and why the evol ved Am stars are now known as p Pup stars. 

iKurtzl (| 19781 ) presented members of the p Pup class that 
also pulsate. A particularly good case-study of a pPup star 
that is als o a high-amplitud e S Set (HADS) star is that of 
HD 40765 (|Kurtz et alj|l995h . There it was found that spec- 
tral peculiarities still remain in the presence of a peak-to- 
peak surface radial velocity range of 14 km s~ . 

1.2 Models of pulsating Am stars 

The introduction of time-dependent convection in the mod- 
elling of 8 Set star pulsations was a major step towards a 
quantitative c ompa r ison b etween theory and observations. 
iDupret et ail (|2004 l2005h have successfully explained the 
red edge of the 8 Set instability strip by using these models. 
T his approach was succ essfully applied to three 8 Scuti stars 
bv lDupret et~ai1 (|2005h . The S Set stars observed by Kepler 
offer a higher-level test for the theory. In this work we apply 
time-dependent convection models to a 8 Set star observed 
by Kepler for the first time. 

Before we move on to the specific case of KIC 3429637, 
let us discuss the issue of pulsation coincident with metalli- 
cism, i.e. of the pPup stars. In the Am stars it was thought 
that there should be insufficient helium left in the He II par- 
tial ionisation zone for 8 Set pulsations to be driven, but this 
is bec o ming increasi ngly doubted (e.g. ICatanzaro fc Balonal 
I2012D . IKurtj (|l989t ) discussed the theoretical exclusion of 
metallicism and pulsation in detail before announcing the 
first case of an extreme classical Am star showing 8 Set 
pulsations. It appears that pulsation in chemically pecu- 
liar stars with pulsation velocities of several hundred ms" 1 
can be so laminar that turbulence on the order of cms" 1 
is not generated. If turbulence were generated, the diffu- 
sion process would be come ineffective and the stars would 
become homogenised. iTurcotte et alj (|2000t ) addressed the 
issue, noting that it is velocity gradients that generate tur- 
bulence, and not the speed of displacement itself - fast but 
uniform displacement will not become turbulent. Their mod- 
els, the New Montreal Models, which take diffusion of el- 
ements up to Ni into account, find stable stars with He 
abundances as low as 0.114 at 750 Myr. Specifically, He 
is still substantially present in the Hell driving region in 
these models, and the iron-peak-element opacity bump is 
increased relative to standard models and contributes to ex- 
citation of longer-period modes. As the stars evolve, they 
begin to pulsate naturally (i.e. without invokin g dredge- up 
or hypothetical mass- loss (|Turcotte et al.ll200Ct l'l because of 
the way that radius expansion causes the Hell ionisation 
zone to shift i nward in mass fraction whe re more residual 
He is present (|Cox. Hodson fc Kind [T979I ). The evolution 



lead s to the increase in period of observed pulsation modes, 
and ITurcotte et al.l concluded that generally the effect of 
diffusion is to stabilise against p-mode pulsations but excite 
g-modes. No extensive observational study covering a large 
sample of stars has been carried out to investigate pulsation 
frequency ranges in light of abundance anomalies. 

In spite of all the strengths of diffusion theory (see 
e.g. lKurtj|2000l ). the diffusion models still fail to repro- 
duce the rece nt observations of th e commonality of pulsat- 
ing Am stars (|Smallev et al.l [20 111 ) satisfactorily, though we 
note that the disparity is not as great as it seems because 
24 per cent of their "pulsating Am" sample were actually 
"Fm 8 Del" stars, and that pulsations are expected to occur 
naturally in these evolved stars according to modern diffu- 
sion mo dels. We refer the reader to the work of lBalona et alj 
|2011bl ) for further discussion and observations on the edges 
of the instability strip for pulsating Am stars. 

1.3 The Kepler Space Mission 

The Kepler Space Mission features an array of 42 CCDs cov- 
ering a 115 deg 2 area of the sky between the constellations of 
Cygnus and Lyra. It captures white-light photometric data 
from ~150 000 stars, with the objective of detec ting transit- 
ing p lanets that orbit within the habitable zone (|Koch et alj 
l20ld ). The attained precision is around the micromagnitude 
level, and, with a duty cycle > 92 per cent, the mission offers 
vast advances in astrophysics through the study of stellar 
pulsations. With photometry alone, transiting studies only 
yield the ratio of the planet and host-star radii, but astero- 
seismology can determine stellar radi i to within ~ 1 per cent 
m some cases |Gilliland et alj|20ld ). With that in mind, 
about 1 per cent of observations are allocated to asteroseis- 
mology. 

Kepler data are available in two cadences: long-cadence 
(LC) with effective 29.4-min integrations; and short-cadence 
(SC) with 58.9-s integrations. The near-continuous obser- 
vations are interrupted every 32 d (one 'month') for data 
downlinking, and every three months corresponds to one 
quarter of Kepler's 372. 5-d heliocentric, Earth-trailing or- 
bit. LC data are therefore available in quarters (denoted 
Qn), and SC data are further subdivided into months (de- 
noted Qn.m). Each quarter the satellite performs a roll to 
keep its solar panels pointed towards the Sun and its radi- 
ator pointed towards deep space. Stars subsequently fall on 
to a different detector module when the symmetrical focal 
plane rotates. The unfortunate failure of module 3 early in 
the mission means that a small fraction (~ of the field 
of view is only observable for three quarters of the year. 



2 PHOTOMETRIC OBSERVATIONS 

Kepler data are available in abundance for KIC 3429637, 
which is the 9th brightest in the 8 Set working group of the 
Kepler Asteroseismic Science Consortium (KASC). Long- 
cadence (LC) public data are available from Ql through to 
Q8, and short-cadence (SC) public data are available for 
each month of Q7 and Q8. The LC data, with just over 
30 000 data points covering 670.33 d, have a duty cycle of 
91.4 per cent. The light-curve for the Pre-search Data Con- 
ditioned (PDC) LC data is displayed in Fig.[TJ from which 
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it is immediately clear that the amplitude of the light vari- 
ation is growing with time. We used only the Least-Squares 
version of the PDC pipeline (i.e. PDC-LS). 

There are two obvious potential sources of the appar- 
ent amplitude growth: (1) the amplitude growth is intrin- 
sic to the star, i.e. astrophysical; and (2) the amplitude 
growth is in some way instrumental. PDC data are not cre- 
ated to facilitate asteroseismology, but rather to prepare the 
data for planet searches. As such, the data are not to be 
used for asteroseismology without caution that the pipeline 
might modify stellar variability. To test whether the pipeline 
caused the amplitude growth we compared the PDC data to 
the simple aperture photometry (SAP) data, for which only 
basic calibration has been performed and which is deemed 
suitable for asteroseismology, provided that instrumental 
trends are corrected for. Since instrumental trends gener- 
ally affect low frequencies only, these are not an issue in the 
periodogram at the frequencies of pulsation shown in the 
Fourier transform of the data in Fig. [2] The amplitudes of 
the pulsations in the PDC data were compared to the uncor- 
rected SAP data, by means of a non-linear least-squares fit 
of the frequencies to each quarter of data u sing the software 
package PERIOD04 (|Lenz fc Bregerl l2004 l. The amplitude 
growth was the same in both data sets. 

The star has two dominant frequencies (/i = 10.337 d _1 
and /2 = 12.472 d -1 ), with amplitudes growing between 1.5 
and 2.5mmag, and a third mode that lies between them 
(/3 = 10.936 d -1 ) with an amplitude of ~0.3mmag - ob- 
servable without any prewhitening. These frequencies are 
typical for 8 Set stars. As an additional check that the am- 
plitude growth is astrophysical and not instrumental, the 
amplitude of that third peak was tracked across multiple 
quarters for comparison. It was found that although the two 
main peaks continued to grow in amplitude across all quar- 
ters, fz showed a continual decrease in amplitude, as seen 
in Fig. [3] This proves that the growth in amplitude of the 
light variation is astrophysical. In addition to these three 
main frequencies, some ~40 statistically significant peaks 
with lower amplitudes exist in the data. We tracked the 
amplitudes of only the next four highest-amplitude modes, 
and have hence tracked all modes for which modelling was 
performed (see §6). The amplitudes of these four modes, 
not shown in Fig. [3] show changes of no more than 30 /imag 
across the 8 quarters and the changes appear to be noise- 
dominated, so the amplitudes of peaks with lower statistical 
significance were not tracked. 

The optimum aperture for a star is defined such that 
the signal-to-noise ratio of the light captured from that star 
is maximised. This doesn't, therefore, capture all light from 
that star, because pixels further from the centre of the star 
contain less of the star's light and therefore relatively higher 
noise levels. It is possible, as a direct consequence of the 
choice of optimum aperture, that slightly more light from 
this star is being captured in each quarter and slightly less 
of a neighbouring star that is contributing towards the total 
light, which would create the effect of amplitude growth of 
all pulsation frequencies in the star's Fourier spectrum. Al- 
though not all modes are growing in amplit ude, a check is 
necessary. The Kepler In put Catalogue (KIC; iLatham et al.l 
l2005l ;l iBrown et al.ll20l"lh lists the contamination for this star 
to be 0.037, meaning that about 3.7 per cent of the light in 
the aperture of this star is coming from one or more neigh- 



bours rather than the target star itself. To rule out the hy- 
pothesis that the amplitude growth is due to changing flux 
fractions of any neighbour, we took the target pixel light 
curves of Ql to Q5 and defined a mask for each one that 
covered every pixel available for the star for that quarter. 
We found that just as for the PDC LS pipeline data, with 
our custom mask the amplitudes of the dominant modes 
grew each quarter and the amplitude of the third mode that 
lies between them continued to decrease. 

In the CCDM (Catalogue of the Compon ents of Double 
and Multiple stars; iDommanget fc Nvslll994h . KIC 3429637 
is listed as a double star with separation 1.58arcsec. The 
primary has V = 7.8 and the secondary has V — 12.7, thus 
the secondary has a spectral type around early K. Having 
a Hipparcos parallax of 3.75 ± 0.58 mas, the stars are well 
separated with an orbital period > 1000 y, so tidal effects 
on the primary's rotation period are unimportant. Without 
radial velocities, we must treat KIC 3429637 as if it were a 
single star. Furthermore, there is no evidence at present to 
rule out chance alignment. 

We wish to examine the possibility that the amplitude 
change is caused by the secondary. According to the CCDM, 
the dimmer star contributes just 1 percent of the light com- 
ing from the pair. For this dimmer star to cause an ampli- 
tude change of ~30 per cent is impossible: one would have 
to remove it thirty times to exact such a change. 



3 SPECTROSCOPIC OBSERVATIONS 

Various methods exist for inspection of a star for chemi- 
cal peculiarity, and definitions of peculiarity in th e liter - 
ature are complex. In his Am star review, IContil (|l970h 
summarised the peculiarities as "deficient Ca (Sc]or over- 
abundant heavier elements, or both." Indeed, IContil (|l965h 
used the line depth ratio of Se ll A424 6 .8 /S rII A4215.5 as 
an indicator of peculiarity and ISmithl (119701 ) went further 
and used another criterion Sell A4320.7/YII A4309.6, be- 
cause of the potential for blending of t he Sell line w ith the 
Fel line (at A4247.3) in fast rotators; ISmithl (| 19701 ) noted 
that a substantial fraction of Am stars may be nearly nor- 
mal in Sc and Ca. Modern abundance surve ys are consis- 
tent w ith long-established characterisations: iGebran et alj 
(I2010D found that: "All Am stars in the Hyades are deficient 
in C and O and overabundant in elements heavier than Fe 
but not all are deficient in calcium and/or scandium." 

Another method for detection of peculiarity is the 
Stromgren Ami index. A-stars with negative Ami indices 
are commonly Am stars, with more negative values indicat- 
ing a more chemicall y peculiar sta r. U sing uvby parameters 
for KI C 3429637 from lOlsenl (Il983h and lHauck fc Mermilliodl 
(|l998h (b-y = 0.189, mi = 0.215, ci = 0.849), noting that 
both sources agree to that precision, t he Ami in d ex for this 
star was calculated using Table 1 of ICrawfordl (|l979h for 
calibration as A mi = —0.028, which indicates an Am/ Am : 
star. Our evaluation is consistent with that of I Abtl (|l984h . 
who classified the star as Am(F2/A9/F3) - a marginal Am 
star given that the metallic-line type and Ca II K-line type 
differ by fewer than five spectral subclasses. We note that al- 
though the Crawford A mi calibration is for main sequence 
stars, it may also be applied to evolved stars. 

Atmospheric parameters for KIC 3429637 were calcu- 
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Figure 1. Upper Panel: the <5 Set pulsations of KIC 3429637 in SC PDC from Q8. Kp denotes Kepler magnitude, and BJD is Baryccntric 
Julian Date. Lower Panel: the unedited light-curve of the LC PDC data from Ql to Q8. The amplitude of the light variation is clearly 
growing with time, as can be seen from the unresolved envelope of the pulsations. 



lated by ICatanzaro et al.l ()201ll ) as part of their spectro- 
scopic survey of potential Kepler asteroseismic targets in the 
5 Set and 7 Dor instability strips. They estimated an equiva- 
lent spectral type of F0III (from fundamental parameters), 
v sin i = 50 ± 5 km s _1 , and used a variety of methods to de- 
termine T e d and log g. The weighted mean of their values for 
T cti an d log g are 7300 ± 200 K and 3.16 ±0.25 (cgs), respec- 
tively. ICatanzaro et al.l also computed masses for this star: 
a canonical value of 3.6lQ g Mq and a non-canonical value 
(that is, with convective overshooting: Aov = 0.2 H p ) of 
3.21q'j Mq. They estimated the luminosity to be 20 ± 6 Lq. 
Additional parameters from the KIC, roughly derived from 
broad-band photometry, are R — 4.1 Rq and [Fe/H] = 0.12. 
It is immediately apparent from the luminosity class and 



log g parameters that the star is a giant, and thus an evolved 
Am, or pPup, star. A co mparison of t he mi value of this 
star with those studied bv lKurtzl (|l976l . his Fig. 3) strongly 
implies this star is of the p Pup type. We return to this in 
Section[4] 



3.1 New spectroscopic observations 

A medium-resolution optical spectrum of the investigated 
object was obtained using the cross-dispersed, Fibre-fed 
Echelle Spectrograph (FIES) on the 2.5-m Nordic Optical 
Telescope, at Roque de los Muchachos, La Palma during the 
run of 2010 Aug 3-5 (proposal: 61-NOT7/10A). FIES offers 
an optical spectrum from 3700 to 7300 A, with a spectral 
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Figure 2. The Fourier transform of the unedited Q5 LC time-string. The three highest-amplitude peaks, in descending amplitude order, 
arc /i, f 2 and / 3 . 
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Figure 3. The growth of the amplitude of /i and /2, and decrease in amplitude of /3 in LC PDC-LS flux as a function of time. Linear 
fits are shown for /i and /3, and an exponential fit is shown for /2. Errors on amplitude are on the order of a few micromagnitudes, and 
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resolution of 46 000 in a single exposure. The signal-to-noise 
of the spectrum is about 100 at 5500 A. 

The spectrum was reduced using standard procedures 
of FIEStool, which consists of bias subtraction, extraction 
of scattered light produced by the optical system, division 
by a normalised flat-field, and wavelength calibration. After 
reduction, the spectrum was normalised to the continuum 
by using SPLAT, the spectra l analysis tool from the Starlink 



by using sFLAT, the spectra l 
project (jDraper et al.lboOSl l . 



3.2 Atmospheric parameters determination 

To perform an abundance analysis, one needs to determine 
an appropriate atmospheric model of the star, which re- 
quires the knowledge of its effective temperature T B g, sur- 
face gravity log g and metallicity. The necessary atmospheric 
mode ls were compute d with the line-blanketed LTE ATLAS9 
code l|Kuruc j Il993ah , which treats line opacity with opac- 
ity distribution functions (ODFs). Th e synthetic sp ectra 
were computed with the SYNTHE code jKuruc z 19931)) over 
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Table 1. The elemental abundances for KIC 3429637 as compared to the solar values of lAsplund et al Columns 4 and 5 display 

the abundances of this star and the Sun, respectively, with their standard deviations. For those elements in column 4 whose abundances 
were determined from only one or two lines, an accurate uncertainty estimate cannot be made so we provide the mean uncertainty of 
our other ab undances (0.12 ) inste ad. Column 6 is the difference between columns 4 and 5. Column 7 contains [M/H] for the normal 
A- stars from [ G ebran et alJ i201Ch for the elements they analysed. In their paper, these abundances were relative t o the solar values 
of iGrc vcssc fc Sauvall limb , but here we correct them to the latest and most precise abundances: lAsplund et al. | j200Sh . Column 8 
contains the difference between the abundances of this star and a normal A-star. The final column is like column 8, but with a metallicity 
correction applied - it represents the difference between KIC 3429637 and a normal star of the same Fe/H (see in-text explanation, 
Scction r3,3H . 



element 


Z 


# of lines 


abundance 


solar value 


[M/H] 


A-stars 


AM 


A M corr 


C 


6 


6 


8.03 ±0.16 


8.43 ± 0.05 


-0.40 ± 0.17 


-0.06 


-0.25 


-0.49 


Na 


11 


1 


6.07 ± 0.12 


6.24 ± 0.04 


-0.17 ± 0.13 


0.29 


-0.37 


-0.04 


Mg 


12 


6 


7.53 ± 0.15 


7.60 ± 0.04 


-0.07 ± 0.16 


0.13 


-0.22 


0.07 


Si 


14 


4 


7.32 ± 0.17 


7.51 ± 0.03 


-0.19 ± 0.17 


0.39 


-0.54 


-0.08 


S 


16 


1 


7.67 ± 0.12 


7.12 ± 0.03 


0.55 ± 0.12 








Ca 


20 


11 


6.34 ± 0.09 


6.34 ± 0.04 


0.00 ± 0.10 


0.04 


-0.02 


-0.03 


Sc 


21 


11 


2.97 ±0.10 


3.15 ± 0.04 


-0.18 ± 0.11 


-0.25 


0.09 


-0.38 


Ti 


22 


29 


4.82 ±0.13 


4.95 ± 0.05 


-0.13 ± 0.14 


0.15 


-0.21 


-0.11 


V 


23 


2 


4.43 ±0.12 


3.93 ± 0.08 


0.50 ±0.14 








Cr 


24 


25 


5.50 ±0.20 


5.64 ±0.04 


-0.14 ±0.20 


0.17 


-0.28 


-0.05 


Mn 


25 


6 


5.24 ±0.14 


5.43 ± 0.05 


-0.19 ± 0.15 


0.13 


-0.36 


-0.16 


Fc 


26 


80 


7.35 ±0.11 


7.50 ±0.04 


-0.15 ± 0.12 


0.19 


-0.34 


0.00 


Ni 


28 


17 


6.43 ±0.12 


6.22 ± 0.04 


0.21 ± 0.13 


0.36 


-0.12 


0.39 


Cu 


29 


2 


3.92 ±0.12 


4.19 ± 0.04 


-0.27 ±0.13 








Zn 


30 


2 


4.37 ±0.12 


4.56 ± 0.05 


-0.19 ± 0.13 








Sr 


38 


2 


3.97 ±0.12 


2.87 ±0.07 


1.10 ±0.14 


0.39 


0.81 


1.41 


Y 


39 


3 


2.94 ±0.01 


2.21 ± 0.05 


0.73 ± 0.05 


0.57 


0.19 


1.01 


Zr 


10 


3 


2.98 ±0.03 


2.58 ± 0.04 


0.40 ±0.05 


0.62 


-0.20 


0.10 


Ba 


56 


3 


3.08 ±0.13 


2.18 ± 0.09 


0.90 ±0.16 








La. 


57 


1 


0.92 ±0.12 


1.10 ±0.04 


-0.18 ± 0.13 









the wavelength ranges 4000-5850 and 6100-6800A, excluding 
those segments contaminated with telluric lines. Both codes, 
ATLAS9 and SY NTHE were ported under GNU Linux by 
ISbordond (|2005l ) and are available onlin^E The stellar line 
identification and the abundan ce analysis were perform ed 
on the basis of the line list from ICastelli fc Hubrid (|2004l ffl 
We derived T c ff for KIC 3429637 using the sensitivity of 
hydrogen lin e wings to temperature, following the m ethod 
proposed by Ivan 't Veer-Menneret fc Megessierl (|l996T ). The 
effective temperature was estimated by computing the AT- 
LAS9 model atmosphere which gives the best match between 
the observed H5, H7, H/3 and Ha line profiles and those com- 
puted with SYNTHE. The same value of T e s was obtained 
from the analysis of iron lines. In this method effective tem- 
perature, surface gravity and microturbulence (£) are deter- 
mined by the comparison of the abundances obtained from 
neutral and ionised iron lines. The analysis is based on Fe 
lines because they are the most numerous in the stellar spec- 
trum. In general, we require that the abundances measured 
from Fel and Fell lines yield the same result. The absorp- 
tion lines of Fel depend mainly on T e ff, £ and metallicity, 
and are practically independent of logy, whereas the Fell 
lines are mostly sensitive to log g. First, we adjust £ until 
we see no correlation between iron abundances and line in- 
tensity for the Fel lines. Second, T c g is changed until we see 
no trend in the abundance versus excitation potential of the 
atomic level causing the Fel lines. £ and T e ff are not inde- 

1 wwwuser.oat.ts.astro.it/atmos/ 

2 http:/ /wwwuser. oat. ts.astro.it/castclli/grids. html 



pendent of each other. Then log g is obtained by fitting the 
Fell and Fel lines and by requiring the same abundances 
from both neutral and ionised lines. 

The abundances of chemical elements were determined 
by the spectrum synthesi s method. Our analysis follows the 
methodology presented in lNiemczura. Morel fc Aertsl (|2009h 
and relies on an efficient spectral synthesis based on a least- 
squares optimisation algorithm. This method allows for the 
simultaneous determination of various parameters involved 
with stellar spectra and consists of the minimisation of the 
deviation between the theoretical flux distribution and the 
observed normalised one. The synthetic spectrum depends 
on the stellar parameters, such as T a a, logg, £, vsini and 
the relative abundances of the elements. Some of these pa- 
rameters have similar influence on stellar spectra and have 
to be known before the determination of the abundances of 
chemical elements. T e g , log g and £ are the input parameters. 
All other aforementioned parameters can be determined si- 
multaneously because they produce detectable and differ- 
ent spectral signatures. The v sin i values are determined by 
comparing the shapes of observed metal line p rofiles with 
the computed profiles, as shown bv lGravl(|2005h . 

For the chemical abundance analysis we used short seg- 
ments of the spectrum to isolate individual spectral features, 
and in a few cases, some blended features. In the case of 
blended lines, more than one chemical element can influ- 
ence the line profile. Only the most important elements for 
characterising Am star peculiarities were considered. Every 
selected part of the spectrum was analysed by the spectrum 
synthesis method described above. We iteratively adjusted 
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Table 2. Fundamental atmospheric parameters determined from 
the FIES spectrum. 



T e ff 




5 


v sin i 


K 


(cgs) 


kms 


kms 


7300 ± 100 


3.0 ±0.1 


4.0 ±0.5 


51 ±1 



radial velocity, vsini, and chemical abundances until the 
determined parameters remained the same within 2 per cent 
for three consecutive iterations, thus determining the clos- 
est match between the calculated and observed spectrum. 
Finally, we determined the average values of radial velocity, 
of v sin i, and abundances of 20 elements with error estimates 
(Table[T]). Table[2] summarises our determined atmospheric 
parameters, and we plot our spectrum against comparisons 
in Fig.H 

3.3 Discussion of spectroscopic results 

The line ratios of IContil (| 19651 ') and ISmithl (|l970t ) were in- 
conclusive on the Am: nature of the star, necessitating the 
full abundance analysis (TableQ]). We note normal Ca and 
nearly-normal Fe, but slightly deficient Sc, strongly deficient 
C, and highly overabundant Sr, Y and Ba that would collec- 
tively confirm t hat the star is chemi cally peculiar. However, 
observations of iGebran et al.l (|2010l ) suggest that even nor- 
mal A-stars are overabundant in Sr, Y and Zr with respect 
to solar values. This effect of overabundant heavy elements 
can have multiple origins: 

1. The study of IGebran et al.l focussed on the Hyades 
cluster, and those cluster stars could be systematically en- 
riched in heavier elements. The abundances of nearly all 
heavier elements are kn own to be correlated with iron 
(|Hill fe Landstreetj|l993l ). The of the Hyad es clus - 
ter is not well agreed upon. iBoesgaard fe Friell (|l990h 
determined a super-solar metal l icity ( [-^-] ) = 0.127 ± 
0.022 dex, vet lVarenne fc Monierl (| 19991 ) found a much lower 
value of —0.0 5 ± 0.03 dex when investigating 29 F-dwarfs. 
IGebran et alj themselves find = 0.19 dex above solar . 
A recent independent calculation (|Carrera fe Pancinoll201lT ) 
based on three K-giants yields ([§]) = +0.11 ± 0.01 dex 
and is close to the mean of previous studies collated by those 
authors. We may therefore safely conclude at least a slightly 
super-solar Fe abundance and that correlation with Fe is a 
contributor. 

2. Using the Sun as a sta ndard star has the we akness 
of assuming the Sun is typical. iMelendez et al.l (|2009l ) found 
the Sun to be depleted in refractory elements by ~20 per cent 
relative to volatile elements compared to solar twins. Specif- 
ically, the lighter, volatile elements (e.g. C, N and O) are 
enriched by ~0.05 dex and the heavier, refractory elements 
depleted by ~0.03 dex. Their findings were confirm ed in an- 
other study ijRamfrez. Melendez fe Asplundll2009l ). and the 
hypothesis that planetary form ation is the cause is indepen- 
dently theoretically supported (|Chambersll2010l ). The pecu- 
liar solar abundance therefore also contributes to the ap- 
parent overabundance of heavy elements in normal A-stars, 
and we may assume that a specific physical effect in A-star 
photospheres is not at play. 

3. One must also consider that the abundances of heavier 



elements are derived from fewer lines, and are therefore at 
greater risk of having erroneous values caused by blends. 

Nevertheless, for those elements that were studied 
by IGebran et al.l we also show the difference between 
KIC 3429637 and normal A-star abundances, relative to the 
Sun in T ableg] Mn, Ni, Sr and B a are str ongly correlated 
with Fe IGebran, Monier fe Richard! |2008|V and Y and Zr 
increase rapidly with Fe/H (jGebran et aljboid ). thus the 
chemical peculiarities of KIC 3429637 would be even more 
pronounced were it not for its low iron abundance. There- 
fore in the final column of Table[T]we tentatively provide the 
abundance differences as if KIC 3429637 had the same Fe/H 
as the mean value of the 'normal' A-stars to which we are 
comparing, namely [Fe/H] = 0.19. The correction is calcu- 
lat ed from the data ta bles provided in the online material 
of IGebran et~all |2010l ). Since A-stars are known to show 
large star-to-star variations in abundan ces compared to F- 
stars (|Gebran. Monier fe Rich ard 2008), the correlation be- 
tween [M/H] and [Fe/H] is not always tight, hence the 'ten- 
tative' provision. Importantly, the elements for which the 
correlation is tighter are those that are the most peculiar. 
iFossati et alj (|2008bl ) found that there is no temperature ef- 
fect on abundance patterns, thus we need not worry about 
any temperature difference between KIC 3429637 and the 
mean of the A-stars against which we are comparing. 

We did not take non-LTE effects into account in our 
calculations, yet these effects can be important. Given the 
importance of calcium in Am star spectra, we considered 
non-LTE corrections for this element , but they are small 
|Mashonkina. Korn fe Przvbillal [20071 ). For iron, non-LTE 
corrections are significan t (0.2 — 0.3 dex) for F el, but not 
Fell, and are positive ijRentzsch-Holml Il996l ). The cor- 
rections would bring our Fe abundance closer to the so- 
lar value. Conversely , corrections for carbon are negative 
llRentzsch- Holmll996l ) and would further exaggerate the car- 
bon deficiency seen in this star and expected in Am stars. 

All three of our T e ff, logg and vsini determinations 
agree extr emely well with t h e weig hted means of those pre- 
sented in ICatanzaro et al.l (|201ll), and £ is slightly high 
for a star of this temperature |Takeda et alj 120081). We 
note, however, that in Am stars £ i s expected to be 
slightly higher than in norm al stars (|Gebran fe Monierl 
120071 : ICouprv fe Burkhartlll992l ) . The latter authors pointed 
out that the £ distribution reaches a maximum around 
A5, where the angular rotational velocity distribution also 
reaches a maximum. This could implicate rotational insta- 
bilities in generating microturbulence, but if that were true 
Am stars would not be observed to have higher £ values. 
IPace et al.l (|2006l ) defined the relation £ = -4.71og(T cff ) + 
20.9 kms -1 for horizontal branch stars, in which radiative 
levitation is also important - this would give an expected 
£ of 2.7kms _1 , but the relation is not appropriate for the 
whole main sequence, and might even be used as a diagnostic 
to distinguish horizontal branch stars from main sequence 
stars. 

The low log g value indicates this star is evolved. In- 
terestingly, the value of v sin i is high for an evolved Am 
(i.e. pPup) star - Am Main-Sequence stars are already slow 
rotators, with the Am/Ap star unimodal rotation velocity 
distribution peakin g at ~60kms~ 1 afte r statistical correc- 
tion for inclination (|Abt fe Morrel]H"l995l ). During evolution, 
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Figure 4. Spectrum of KIC 3429637 (centre) with two comparisons: the FOIII star HD7312, and the F5III star HD 65925, though 
neither is an MK-standard. ICatanzaro et al" . (2011) listed KIC 3429637 as F0 III m. The hydrogen line type is a good match to 
~F0/A9. Despite being of the Am type, notice that the Call K-line is no weaker than the F0 comparison (cf. Tablcfl] in 
which calcium is completely normal). The metal lines are much stronger, though, and are even stronger than the F5 compar- 
ison. In particular, the SrII A4077 line is much deeper, owing to the Am nature. The comparison spectra were obtained at 
'http://www.cso.org/sci/observing/tools/uvespop/ficld_stars_uptonow.html', normalised at 3910. OA and vertically separated by 1 unit. 



stars transfer angular momentum inward through radius ex- 
pansion, slowing the surface rotation further. Hence by ac- 
cepting the star is evolved, the high usini argues against 
the type of close binary system common among Am star^f). 
Instead, the hypothesis of a single A star with a sub-average 
rotation velocity is favoured for this case, in which this star 
was rotating slowly enough to develop some peculiarities. 

It is well established that a large fraction of Am stars 
are found in binaries and that tidal braking slows the stel- 
lar rotatio n, allowing diffusion to occur. With the exclu- 
sion of Sc |Fossati et al.|[2008al ). pecu liarities are more ex - 
treme in more slowly rotating stars (jTakeda et al.|[200ct ). 
The observed moderate v sin i value is indeed compatible 
with the fact that KIC 3429637 is not an extreme Am star, 
as the moderate abundance anomalies imply. Notwithstand- 
ing, pPup stars do not have extreme abundances , regardless 
of v sin i. However . I St at e v&, Iliev fc Budail (|2012h found that 
peculiarity is also slightly correlated with orbital eccentric- 
ity: more eccentric orbits produce more anomalous abun- 
dances, but we do not expect this to have any significance 
for wide binarie s. Attempts to use the freq uency modula- 
tion technique of lShibahashi fc Ku rtz (2012) to infer orbital 
parameters were unsuccessful - orbital sidelobes were not 
present. This indicates that either the star is in a binary 



3 Tidal braking is most significant when the orbital period de- 
creases belo w 10 d - the tidal braking mechanism proposed by 
IZahnl lll975l ) denotes that braking timescales depend on the stel- 
lar separation to the power 8.5. This is confirmed observationally 
by the fact that S B2 systems are circular for periods below ~8 d 
iDebernardill2000l) . 



system with an orbital period longer than the dataset, i.e. 
forb > 700 d, or that the star is not binary at all. 



4 ASTROPHYSICAL CAUSE OF AMPLITUDE 
GROWTH 

The different amplitude growth rates of each mode, visi- 
ble in Fig. [3] can be used to determine a time at which the 
mode amplitudes were too low to be detectable from the 
ground, assuming the growth rates are well behaved. A lin- 
ear growth rate gave the best fit to fi, but an exponential 
growth rate was more appropriate to y_, which is growing 
much more rapidly than the lower frequency mode. Taking 
the limit of ground-based detection as 1 mmag in the ampli- 
tude spectrum, /1 would only have become detectable ~2.6 y 
ago, and _/_ ~11 y ago. We might c onclud e that these modes 
would have been undetectable in Il984l at the time when 
lAbtl spectroscopically observed the star to be metallic lined, 
although no photomet ric inve s tigati on into variability was 
carried out. However, iBregerl (|200Ch observed the evolved 
8 Set star 4 CVn over three decades and found mode am- 
plitudes to be unpredictable and highly variable. One high- 
amplitude mode in that star disappeared altogether and re- 
emerged with a random phase. With so few studies of am- 
plitude growth rates, and the unpredictable behaviour seen 
in 4Cvn, we cannot predict how KIC3429637 will behave 
or how it behaved in the past. This star is most certainly 
worthy of continued observation to monitor the amplitude 
changes. 

In any case, the current rates of growth indicate the 
star has either: (a) not been pulsating for very long, and 
has just crossed some instability threshold; (b) undergoes 
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large amplitude variations over long periods of time, possi- 
bly in some cyclic nature; or (c) none of the above. We do not 
have enough data to pursue hypothesis (b), in that we have 
not observed a full cycle of variation in the light curve, but 
since the amplitude growth is not the same for each mode, 
nor does it appear to be sinusoidal, the cyclical amplitude 
variation hypothesis is unsupported by the present evidence. 
If we assume this star to be like 4 CVn, then continued ob- 
servations are required to analyse the changing amplitudes. 
In the meantime, let us examine hypothesis (a) from a the- 
oretical point of view. 

Even if the pulsations are not especially laminar, the 
destruction of the chemical abundance anomalies may take 
time, such as in the case of oLeo A, which is an evolved Am 
star whose abundance anomalies have not yet been erased 
(|Michaud. Richer fc Richardll2005l ). KIC 3429637 could well 
have been a classical Am star on the main sequence that is 
now at or approaching a stage of rapid evolution in which 
its chemical peculiarities are being erased and its pulsa- 
tion amplitudes are growing. Temporal evolution of abun- 
dance anomalies is another under-studied area of Am star re- 
search; KIC 342 9637 would make a perfect case -study. Exist- 
ing studies (e.g. iGebran et aHl20lA lAbtlll979l) have mostly 
focussed on multiple co-eval stars in one cluster and con- 
trasting clusters of different ages. In these cases the scatter 
in abundances of given elements is large enough that no real 
difference of any element with age is prominent. Focussing 
on single, regularly-observed targets such as KIC 3429637 is 
thus potentially astrophysically rewarding. 



5 NOISE IN THE RESIDUALS 

Spectroscopic data show KIC 3429637 to be cool for 
a 5 Set pulsator, and progressively cooler A-type stars 
have increasingly deep su rface convection zones (see 
iKallinger fc Matthews! |2010| for a discussion). We used the 
Q7.1 SC data to look for such a signature of granulation 
in the form of a rapid change in noise with frequency. We 
present the results in Fig. [5] After extensive investigation, 
the rapid decrease in noise seen between 20 and 30 d -1 
was determined not to arise from granulation, but rather 
the noise is injected by the Kepler PDC- pipeline for 
this particular Kepler quarter. It seems that approximately 
15 per cent of stars in the <5Sct instability strip are afflicted 
with the same problem, which affects both long- and short- 
cadence data, but not necessarily all quarters for a given 
star. A more detailed description will be given in a future 
publication (Murphy, in preparation). 

The Q8 data for this star contain no evidence for the 
"granulation signature", implicating noise injection by the 
pipeline over an astrophysical origin. The drop in noise at 
~24d _1 in SC data is characteristic for this noise injection. 
The Nyquist frequency of LC data prevents the drop being 
seen in that cadence, but the noise is still present. The new 
PDC-MAP (Maximum A Posteriori) data will fix this for 

4 "LS" describes the Least Squares algorithm used by the 
pipeli n e. We refer t he re ader to the papers of IStumpe et alj 
i2012l ), [Smith et al.l [I2012T ) and the Kepler Data Chara cteristics 
Handb ook for more information on the pipeline, and to iMurphvl 
for its effect on asteroseismic analyses. 



the LC data when the old data are reprocessed, but for now 
there is no implementation of MAP for SC. 



6 STRUCTURE AND OSCILLATION MODELS 

In our theoretical computa tions we used the Co de Liegeois 
d'Evolution stellaire CLES (|Scufiaire et al.ll2008T l. where the 
input physics is as follo ws: the equation of state is the su- 
perio r 'CEFF' (see, e.g. IChristensen-Dalsgaard fc DaepperJ 
1992h; the opac i ty ta bles are OPAL opacities from 
Iglesias fc Rogers! (|l996t ) and lAlexander fc Ferguson! (|l994T l 
for high and low temperatures, respecti vely. The relative 
mixtu re of chemical composition from iGrevesse fc Noelsl 
|l993l ) was used, while t he convection was treated us- 
ing Mixing-Length-Theory (|Bohm- Vitensdl 19581 . MLT). Our 
models include neither rotation nor diffusion. For our non- 
adiabatic, n on-radia l oscil lation calculations we used the 
MAD code (|Dupretl [200 j ). where Gabriel's treat ment of 
time- dependent convection has been implemented ll Gabriel! 
119961) . 

We consequently produced a grid, where £ = — 3 eigen- 
frequency spectra are computed for each of the models. For 
the selection of the best model we used the max i mum like- 
lihood method as described by iGrigahcene et al.l (|2012T ). 

6.1 Results and discussion 

Using Kepler data, lUvtterhoeven et al.l (|201lf ) showed how 
A-F pulsators can have such complex frequency spectra that 
they are challenging targets for any asteroseismic study. We 
consequently need as many constraints as possible for any 
theoretical modelling. Although the number of detected fre- 
quencies is reduced in KIC 3429637 compared to what is 
usually detected in this type of star, its modelling is no less 
complicated. 

We performed careful and detailed modelling, varying 
all the parameters (mass, metallicity, mixing-length param- 
eter omlt, etc.) to obtain the combination giving the maxi- 
mum value of the likelihood function. We explored different 
possibilities, taking into account the detected frequencies 
and their associated amplitudes, and the position of the star 
in the HR diagram. The main requirements on our models 
were: 1) fitting the maximum number of the observed fre- 
quencies, 2) modes should be unstable, and 3) the global 
parameters should be inside the observational photometric 
error box. In Table[3]we show the main properties of the best 
obtained model. All the parameters are within the spectro- 
scopic uncertainties. We show the star's position on the HR 
diagram and its evolutionary stage in Fig. [6] 

Some important results come out from our comparison 
between theoretical time-dependent convection models and 
observations. First, we managed to obtain a model where 
seven of the eight highest-amplitude frequencies are excited. 
The lowest frequency mode (fs=0. 39448 d _1 ) is predicted to 
be stable (i.e. not excited). This mode might not be caused 
by oscillation - a t these low fre quencies, instrumental ori- 
gins are common (|Murphvll2012l . his Fig. 3). The origin does 
not appear to be rotational or equal to the difference in 
frequency between any two other modes in Tabled The al- 
ternative to instrumental origin is that it is indeed a pul- 
sation frequency, in which case the conclusion is that some 
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Figure 5. There is a rapid decrease in noise level at around 24 d . Excess noise is injected by the PDC-LS pipeline at frequencies below 
this limit. There are no peaks remaining below that frequency with signal-to-noise > 4.0, and no peaks have been removed at all above 
24 d — 1 , hence the effect has not been caused by preferential prcwhitening. Note that here we have prewhitened many more frequencies 
than were modelled in order to demonstrate the effect more clearly. 



M T cff Log(L/L Q ) Log g R Age X Z £*mlt a 0v 

(M ) (K) (R Q ) (y) 



2.178 7452 1.570 3.648 3.666 8.982E+08 0.7360 0.0149 2.0 0.2 



Table 3. Properties of the best model. 



1.8 




1.2 1 1 

4.02 4 3.98 3.96 3.94 3.92 3.9 3.88 3.86 3.84 

Log(T eff ) 

Figure 6. Position of the best model on its evolutionary track in 
the HR diagram, showing the evolutionary stage of the star. As 
this is a model, there are no error bars. 



physics is missing from the model. We did not model all 40+ 
statistically significant frequencies due to the overwhelming 
complexity involved in doing so. 

Second, we found that the value of the mixing-length 
parameter giving the maximum value for the likelihood 
function is qmlt=2, i. e. different from the solar value in 
l|Grigahcene et al1l2012t ). It is known that lvmlt should vary 
from star-to-star and also with temperature as the internal 
structure changes from a solar-like one to those of A- and 
F-stars. As this star is the first S Set star to be treated with 
TDC models, coverage of lvmlt for T c a > 6000 K is currently 
unavaila ble, but for cooler stars we r efer the reader to the 
works of iGrigahcene et al.l (|2012t ) and lTrampedach fc Steinl 
|201lf ) for deeper discussion. 

A third remark regards the fitting of the frequency val- 
ues. In Table[4] we list the theoretical frequencies together 
with the observed frequencies and mode identifications. Four 
of the observed frequencies are very well fitted by the the- 
oretical modes - the difference is less than 0.03 d _1 (0.35 
/iHz). For the other three the difference is less than 0.2 d . 
However, an ambiguity arises for some modes in that there 
are multiple mode identification possibilities of similar likeli- 
hood. This indicates that ground-based multi-colour photo- 
metric or spectroscopic observation campaigns are necessary 
to complete the mode identification. 

Finally, while the asteroseismic effective temperature 
agrees with the spectroscopic one to within less than 2 a, 
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Table 4. Identification of oscillating modes. Negative n values 
denote mixed-modes with g-mode-like character. An additional 
observed frequency, fg=0. 39448 d _1 , was modelled but was not 
found to be excited. Formal least-squares frequency precision is 
3 X lO" 6 d- 1 for fi and f 2 , and 2 X 10~ 8 d" 1 for / 3 . Since our 
models do not include rotation to lift the degeneracy, all modes 
presented are m = 0. 

the log g value from the models (3.65) does not agree partic- 
ularly well with that from sp ectroscopy (3.0 ±0.1). S imilar 
model behaviour was seen in iGrieahcene et al.l (j2012T l. 

6.2 Conclusions from modelling 

We have compared the observed frequencies of KIC 3429637 
to state-of-the-art <5Sct models. The results of this com- 
parison are that: 1) all the observed modes (those listed 
in TableQ are predicted to be excited in our model ex- 
cept the lowest frequency, 2) four out of seven frequen- 
cies are fitted to within 0.26 per cent, and 3) the varying 
amplitudes between quarters might be explained by mode 
interaction - the modelled frequencies are all rather close 
and are perhaps subject to r e sonan ces as described in, e.g. 
iDziembowski fc Krolikowskal |l985|). We identify many po- 
tential mixed modes, and it is noteworthy that g-modes that 
aren't even visible in the Fourier transform have the poten- 
tial t o affect p-mode frequencies jBuchler. Goupil fc Hansen! 
1 19971 ). 



7 CONCLUSIONS 

The high-precision Kepler data indicate two dominant 
modes of growing amplitude, with different growth rates, 
and a third mode of intermediate amplitude whose ampli- 
tude decreases. The mode amplitude growth appears to be 
intrinsic to the star. 

Our spectroscopic data confirm literature classifications 
that this star is of the marginal Am type. The heavy ele- 
ments are overabundant with respect to solar values, and 
the light elements are underabundant - as one expects for 
Am stars. Notably, calcium is normal. 

Mode amplitudes in 5 Set stars can change by large 



amounts on short timescales, with changing growth rates 
(cf. 4 CVn). The amplitude growth in this star appears much 
less chaotic than for 4 CVn, and we have proposed evolution 
as the dominant cause, in that pulsation amplitudes are ex- 
pected to increase naturally in evolving Am stars as the 
He II convection zone deepens and picks up residual helium. 
We cannot definitively rule out that mode interaction is the 
cause, where energy is being transferred from one or more 
modes into those whose amplitudes are growing, except to 
say that we see no modes decreasing in amplitude by com- 
parable amounts to those by which the dominant modes are 
growing. As such, this explanation seems viable only if one 
or more unseen modes - perhaps g-modes or high-degree 
p-modes - are transferring the energy. 

For the first time, time-dependent convection has been 
used in models of a 5 Set star observed with Kepler. Despite 
the narrow observational constraints, the modelling was dif- 
ficult. We determine a 2.18-Mq, 7450-K model with Z = 
0.0149 to be the best solution (cf. Table[3|, with the star be- 
ing 94 per cent of the way through its main sequence lifetime 
(see Fig.[5|, and we match four of seven mode frequencies 
to within 0.26 per cent (Table|4|. Multi-colour photometry 
and/or time-series spectroscopy would facilitate mode iden- 
tification, and this target is most worthy of continual follow- 
up in Kepler's long-cadence mode, not least because it offers 
a chance to witness stellar evolution on human time-scales. 
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